Charge qubit dynamics in a double quantum dot coupled to phonons 
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The dynamics of cliarge qubit in a double quantum dot coupled to phonons is investigated the- 
oretically in terms of a perturbation treatment based on a unitary transformation. The dynamical 
tunneling current is obtained explicitly. The result is compared with the standard perturbation the- 
ory at Born-Markov approximation. The decoherence induced by acoustic phonons is analyzed at 
length. It is shown that the contribution from deformation potential coupling is comparable to that 
from piezoelectric coupling in small dot size and large tunneling rate case. A possible decoupling 
mechanism is predicted. 
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I. INTRODUCTION 

Computers based on quantum mechanics is proven, by 
computational theory, more efficient at some specific cal- 
culations than those based on classical physicsiSi^ The 
first step to build a quantum computer is the realiza- 
tion of its building block, quantum bit (qubit). Within 
the last decade, various schemes have been proposed and 
many of them have even been realized Among 
them, the electrically controlled charge qubit in a semi- 
conductor double quantum dot has the potential advan- 
tages of being arbitrarily scalable to large system and 
compatible with the current microelectronics technology. 
Besides, the double-dot system is also extremely useful 
in basic physics as it enables us to investigate the deco- 
herence and dissipation of a small quantum system inter- 
acting with its environment. 

Various designs of double-dot qubits have been studied 
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Recently, Hayashi et al. success- 



fully realized coherent manipulation of electronic states 
in a double-dot system implemented in a GaAs/AlGaAs 
heterostructure containing a two dimensional electron 
gas. The damped Rabi oscillation was observed in the 
time domain in their experiment and an experienced 
formula was presented to fit the experimental data^ 
Fujisawa et al. explained the transport in this system 
through density matrix simulations.^"* But the phonon 
effect was not included. However, the phonon induced 
decoherence is significant according to the analysis of 
Ref. [13]. Considering the interaction with the piezo- 
electric acoustic phonons, Brandes et al. investigated in 
detail the single electron tunnelling in a double quan- 
tum ciot>i^*iSiiLi2iiS Fedichkin and Fedorov also stud- 
ied the error rate of the charge qubit coupled to an 
acoustic phonon reservoir ^2i2i But no dynamical tun- 
neling current is presented by all these work. In a latest 
review paper,^^ Brandes derived the dynamical tunnel- 
ing current in the weak electron-phonon dissipation limit 
through Born-Markov approximation (PER). But the re- 
sult is rather complicated. 

In this work, we study the damped Rabi oscillation ob- 
served in Ref. [13]. The quantum dynamics of a single 
electron tunneling in the double-dot system is investigate 



without applying the Born-Markov approximation to the 
electron-phonon interaction. A simple explicit expression 
of dynamical tunneling current is presented through per- 
turbation treatment based on a unitary transformation. 
The phonon (both the deformation potential and piezo- 
electric contribution are included) induced decoherence 
is investigated at length and possible decoupling mecha- 
nism is presented. 

This paper is organized as follows. In Sec. II we in- 
troduce the model Hamiltonian for the charge qubit and 
solve it in what follows. A comparison with PER ap- 
proach is given in Sec. HI. In Sec. IV, we derive the 
spectral density of a double quantum dot and then an- 
alyze the phonon induced decoherence. The conclusions 
is given in the Sec. V at last. 



II. MODEL AND THEORY 

In this section, we shall introduce the model Hamilton- 
ain for the double quantum dot and then develop general 
theory leading to explicit expression for the dynamical 
current. We will focus on the phonon effect on qubit 
dynamics. 



A. Model Hamiltonian 

The double quantum dot consists of left and right dots 
connected through an interdot tunneling barrier. Due to 
Coulomb blockade, only one excess electron is allowed 
to occupy the left or right dot, which defines two basis 
vectors of \L) and \R) (with the energy level El and er 
respectively) in the Hilbert space. The energy difference 
between these two states e = El—er can be controlled by 
the source-drain voltage Vscf^ Considering the coupling 
to its environment, the double dot can be described by 
the Hamiltonian: 
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H„ + Hp 



Here the qubit Hamiltonian reads(/i = 1) 
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where is the interdot tunnehng, Ux and are Pauh 
matrix with cr^ = \L){L\ - \R){R\ and g.^ \L){R\ + 
\R){L\. If the qubit isolates from any other degrees 
of freedom, the excess electron would oscillate coher- 
ently between two dots with the Rabi frequency A = 
\J -\- 4T^. However, the qubit has to couple to its en- 
vironment (phonons and electron reservoir in leads) in 
practice. Hp and Hf,p stand for the phonon reservoir and 
its coupling to charge qubit respectively. They can be 
written as follows: 



Hp = ^i^c^.\b\^xK^^ (3) 
Hep - \<r.Y.^M^A,x + K..>-Kxl (4) 



where 6^ _^ (^q,A) and uj^^.x are the creation (annihilation) 
operators and energy of the phonons with the wave vec- 
tor q and polarization A, Mq > is the electron-phonon 
coupling constant. The effects of the phonon bath are 
fully described by a spectral density 



E 

q.A 



M. 



(5) 



Hr and H^r in the Hamiltonian H stand for the elec- 
tron reservoir in leads and its coupling to charge qubit 
respectively. 

In experiment, a pulse technique is used to switch the 
the Vsd from large bias in the initialization process (an 
excess electron localizes in the left dot) to the zero bias 
in the manipulation process (the double dot is isolated 
from leads, and the excess electron tunnels resonantly 
(i.e., e = 0) back and forth between two dots).^^ Restor- 
ing a large bias voltage Vsd after the pulse time t gives 
the measurement of dynamical elastic tunneling current 
which stands for the probability n{t) of the excess elec- 
tron in the right dot at that exact time. 

Neglecting the higher order tunneling (cotunneling) 
between leads and the dots, the effective Hamiltonian 
in the manipulation process reads: 



Here, for the sake of simplicity, we omit the polarization, 
since it makes no difference in the theory below. When 
it makes difference (in Sec. IV), it will be included again 
and noted out explicitly. This effective Hamiltonian is 
the starting point for our theory. 



B. Theory 

The effective Hamiltonian H^ff is equivalent to the 
spin-boson Hamiltonian in zero bias case. Though it 
seems rather simple, it cannot be solved exactly. Various 
analytical or numerical approaches have been proposed 
to obtain an approximate solution to itiSi^ 



Here, we apply a canonical transformation: H = 
exp(S')-ffe// exp(— 5) , with the generator t^SiSSiSi 

q ^'^q 

Thus we get the Hamiltonian H , and decompose it into 
H' = Hf^ + H[ + H'2, where 

h'o = r/Tef7,+^c.q6tfc^-^^eq(2-eq),(8) 
q q 

h[ = ia,^(l-eq)(A/q6^ + Af*6q)-Kr;r,*a,A(9) 



H2 = Tcax{coshA-T]) +Tci(7y{sinhA-riA), (10) 



where 



(11) 



and 77 is a parameter which will be adjusted to minimize 
perturbation terms {Hi and i?2)- Obviously, Hq can be 
solved exactly. We denote the ground state of Hq as 



Iff> = k2)|{0q}), 

and the lowest excited states as 

|es) = ki)|{Oq}), 

|eq) = k2)|lq), 



(12) 



(13) 
(14) 



where |si) and |s2) are eigenstates of {(TxIsi) ~ |si), 
CT3;|.S2) = — |s2)), |{0q}) stands for the vacuum state for 
phonon, and |lq) means that there is only 1 phonon for 
mode q and no phonon for other modes. Let Hi\g) — 
and {g\H2\g) = 0, we will get i^q and 77 respectively as 
follows: 



Cq 



cjq + 2riTc ' 



(15) 
(16) 



Now we can easily check that {es\Hi\es) = 0, 
(eq|i?;|eq) - 0, (e.|i?;|.g) = 0, (e^lH^lg) = 0, and 
{eq\H[\es) = Vq, where Vq = With these 

relations above, we can now diagonalize the lowest ex- 
cited states of H as 

h' = -vTc\g){g\ + Y,E\E){E\ 

E 

-|-terms with high excited states. (17) 

The experiment in Ref. [13] is performed at lattice tem- 
perature below 20 mK.^^'^"* At such a low temperature, 
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the multiphonon process is weak enough to be neghgible. and the E^s are solutions to the equation 
So we can get the transformation as^SiSLS^ 



where 



|e.) = 5]x(i?)|£;), 

E 
E 



|x(i?)p = [l + ^ 



{E + 77TC - Wq)2 



|yq(^^)l 



[E + r,T, - Wq)2 



(18) 
(19) 
(20) 

(21) 



E + ?]Tc - UJq 



(24) 



The population inversion can be defined as P{t) = 
{ip{t)\az\tlj{t)) , where \ip{t)) is the total wavefunction 
(qubit and reservoir) in Schrodingcr picture, and 



(25) 



(22) Since the qubit is initialized at the state |L), it is rea- 
sonable to choose IV'(O)) = e"-^|L)|{Oq}). Then we can 

(23) obtain 



q 



(26) 



where uj — E + i]Tc. Denoting the real and imaginary 
part of |V^p/(tj ± iO+ — Wq) as R{uj) and T7(^) re- 
spectively, we can get 



Rio;) 



7(w) 



M7]T,yp / 

Jo 



J(^) 



{uj-uj'){uj' +2r]Tc)^ 
7r^|VqpJ(^-^q) 



,(27) 



47r(r7r,)2 



(28) 



where P stands for Cauchy principal value, and the spec- 
tral density J{uj) is defined in Eq. (5). The parameter 
r] determined by Eq. (15) and Eq. (16) can also be ex- 
pressed as 



?7 = exp-/ duj— — . 

Jq 2[Ll! + 2'qTcY 



(29) 



■ 2{uj + 2riTcY 

The contour integral in Eq. (26) can proceed by calcu- 



lating the residue of integrand and the result is 

P{t) = - cos{ujrt) exp(-7i), (30) 
where we have applied the second order approximation^^ 



7 ~ 7(277Te) = i^J(27/Te), 

and LOr is the solution to the equation 

UJ - 2'nTc - R{uj) = 0. 



(31) 



(32) 



Finally, the tunneling electron population (probablity) in 
the right dot at time t is given by 



7i(t) = i [1 + Pit)] = J [1 - cos(c^,i) exp(-7i)] . (33) 



Thus a rather simple expression for the dynamical tun- 
neling is obtained analytically. The damped oscillation 
indicated by this expression agrees with the experiment 
inRef. [13]. 
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FIG. 1: The tunneling electron population in the right dot as 
a function of time for Ohmic dissipation at low temperature 
limit with the coupling constant g — 0.1, 0.2, and 0.5. The 
resonant tunneling rate is fixed as Tc = O.lWc. The solid 
curve and dot curve stand for the results of our approach 
and PER approach respectively. 



FIG. 2: Long time behavior of the tunneling electron 
population in the right dot as a function of time for Ohmic 
dissipation at low temperature limit with the coupling 
constant g—0.2, 0.1, and 0.05. Others are the same as Fig. 1. 



III. COMPARISON WITH PER APPROACH 

As mentioned in the introduction above, Brandes has 
also derived an expression for the dynamical tunneling 
in the double-dot system. In his derivation the Born- 
Markov approximation is made and the electron-phonon 
coupling Hep is treated as a peterbation (PER)i— In this 
section, we give a comparison between the PER and our 
approach. 

The spectral density J{uj) defined in Eq. (5) is the 
only quantity describing the interaction between the sys- 
tem and its environment that enters into the dynamical 
tunneling. To get the final result, we need the knowl- 
edge of this spectral function first. Ordinarily, it can be 
written in a power law form with a cutoff coc, 

J{uj)=gLo^e--'-% (34) 

where 0<s<l,s = l, and s > 1 corresponds to the sub- 
Ohmic, Ohmic, and super-Ohmic spectral, respectively, 
g is a coupling constant. Since the result for n{t) is only 
presented for the Ohmic case in Ref. [22], we also use 
the Ohmic spectral in this section for convenience. We 



should point out that, our approach is not limited by the 
form of J{i-u). However, the Born-Markov approximation 
made in PER approach is only meaningful and defined 
for s>lM 

Fig. 1 presents the time evolution of the tunneling elec- 
tron population (probability) in the right dot with three 
different coupling constants. The result of both our ap- 
proach and PER is given at the temperature T=20 mK, 
at which the experiment is performed^ . The result at 
higher temperature can also be obtained by PER and we 
find that the difference is indistinguishable when the tem- 
perature is below 0.1 K, which testifies the validity of the 
single phonon process assumption in the transformation 
in Eq. (18), (19), and (20). The comparison shows that, 
in weak coupling regime {g < 0.2), the results of two 
approaches agree with each other well; while in strong 
coupling regime {g ~ 0.5), the difference between two ap- 
proaches becomes clear. Although both two approaches 
are based on perturbation theory, the perturbations in 
two approaches are different. The PER approach treats 
the coupling H^p to the phonon reservoir as a perturba- 
tion, which consequently requires the coupling constant 
g to be restricted in the weak coupling regime. In our ap- 
proach, the perturbation theory is exerted to the Hamil- 
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tonian after a canonical transformation. The perturba- 
tion is taken as and H2, which can be minimized by 
the variational parameter rj. So our approach can be eas- 
ily extended to the strong coupling regime (it works well 
for the whole range of < g < 2).^^ 

Then we focus on the long time behavior of the time 
evolution. Since the double dot is symmetric (i.e., un- 
biased case, £ = 0), the population in two dots at long 
time limit should also be symmetric. In other words, 
the asymptotic value of n{t) is expected to be 0.5. But 
in the strong coupling case {g — 0.5) of Fig. 1, it has 
shown that there is a gap between the result of PER at 
long time limit and that of our approach which equals 
0.5. Actually, even in weak coupling regime, this kind 
of deviation still exists but rather small, as can be seen 
from Fig. 2 in an enlarged view for coupling constant 
g = 0.2, 0.1, and 0.05. Since the asymptotic value from 
non-Markovian calculation presented by DiVincenzo and 
Loss also equals 0.5,^* the deviation in PER approach 
must arise from Markov approximation. And it is avoided 
in our approach. 

There are two alternative ways of perturbation theory: 
one is in the coupling H^p (PER) , another is in the inter- 
dot coupling Tc (POL),^^ To get a good approximation, 
one must restrict the PER approach in the weak cou- 
pling regime, as analyzed above. But it does not mean 
there is no requirement for Tc. In Fig. 3, we present 
the dynamical population in the right dot n{t) in weak 
coupling [g = 0.1), but with three different int-dot tun- 
neling rates Tc = 0.05, 0.01, and 0.001 uJc- We find that, 
in strong inter-dot tunneling case (Tc > 0.05 lUc), two 
approaches give consistent results; while in weak tunnel- 
ing case {Tc < 0.01 Wc), the results deviate from each 
other. It is well known that the energy difference for 
the two eigenstates (i.e., the bonding and anti-bonding 
states) of the double dot system is the Rabi sphtting 
A = a/e^ _|_ 42^2 ^ jjgj-g the unbiased case (e = 0), the 
energy difference A is determined only by the inter-dot 
tunneling Tc- When Tc is too small, the energy difference 
A between two states may be smaller than the pertur- 
bation Hep of the PER approach. In such situation, the 
energy states are nearly-degenerate and the perturbation 
theory for a nondegenerate state of the PER approach is 
not suitable. So the PER breaks down in weak inter-dot 
tunneling regime. In our approach, however, the pertur- 
bation terms {Hi and H2) decrease with Tc, thus the per- 
turbation theory always works well for the whole range 
ofTc. 



IV. PHONON INDUCED DECOHERENCE 

One of the central points in quantum physics is the 
loss of coherence of the quantum system. In this section, 
the spectral density for the double-dot system coupling to 
bulk phonons is derived from the microscopic view. Then 
the phonon induced decoherence is analyzed at length. 
And the suppression of decoherence is discussed. 
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FIG. 3: The tunneling electron population in the right dot as 
a function of time for Ohmic dissipation at low temperature 
limit with the tunneling rate Tc = 0.05, 0.01, and O.OOlf^c. 
The coupling constant is fixed as g = 0.1. 



A. Spectral density 

In last section, we use the Ohmic spectral to calculate 
the dynamical tunneling current. Obviously, it is not 
appropriate to describe the properties of the double-dot 
system. Since all relevant information of the environ- 
ment and its coupling to the qubit system is contained 
in the spectral function, it is necessary and important 
to obtain a suitable spectral function to describe a dou- 
ble quantum dot. Here, we only consider the coupling 
to the bulk acoustic phonons, because we are interested 
in the low temperature limit case. The only two types 
of interaction between electrons and acoustic phonons in 
semiconductors are piezoelectric coupling and deforma- 
tion potential coupling. 

We assume the electron wave functions \L) and \R) are 
sharply around the center of dot with Gaussian shape ~ 
exp[— r^/(2^^)], where I is the the dot size. One can show 
that the piezoelectric coupling constant for GaAs (zinc- 
blende structure) is^ 



1 



2pqsxV 



x(e^2e3 -1-^2^6163 -He3eie2)sin(^),(35) 



where p is the density of the crystal, V is normalizing vol- 
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ume, s is the sound velocity in crystal (longitudinal sound 
and transverse sound have different velocities), = Qi/q, 
^ is the polarization vector whose components depend on 
the polarization mode A, M is the piezoconstant, and d 
is the ccntcr-to-center distance between two dots. Here, 
we include the polarization again, since different polar- 
ization modes give different coupling constants. With the 
simple dispersion relation ujq,\ = s\q, one can now calcu- 
late, according to the definition in Eq. (5), the spectral 
function Jp^(w) due to piezoelectric coupling. But the 
expression is rather complicated. To obtain a tractable 
form of the piezoelectric coupling we use the angular av- 
erage following Bruus et aL'^° and Brandes et ali^. Then 
we get 



0.04 



^sm^)e-'/^^.\ (36) 



where LOd = s/d, loi = s/l, and = P jlii^ ps^ with 

P^M^(i^ + li^), 
35 X 35 



(37) 



where the transverse sound velocity is expressed as x 
times the longitudinal sound velocity s. For the de- 
formation potential coupling, the contribution from TA- 
acoustic phonons is small enough to be neglected as com- 
pared with that from LA-acoustic phonons. So the cou- 
pling constant can be written as^fi 



1 



IpqsV 



yi'e-'^ " /"sin(:V^), 



,q_d, 
2 



(38) 



where is the deformation potential. Then we can easily 
get the spectral function due to deformation coupling 



g,,^^il-^sin-)e-'/'-^\ (39) 

UJ UJd 



where ^df = 

With the parameters of GaAs in Ref . [30] , we can esti- 
mate that (7pz « 0.035 and ^df ~ 0.029 (ps)^^. Previous 
work states that the contribution from deformation po- 
tential phonons is small as compared with piezoelectric 
phonons in double-dot system of GaAs material.-? Our 
result also proves it to be true in the weak confinement 
regime (large dot size). But it is not valid when the dot 
size is decreased to the strong confinement regime. Fig. 4 
shows the spectral functions Jp^(w) and J'^^(w) in strong 
confinement regime, with uji — 1 (ps)~^ (i-e., dot size 
I — 5 nm) and w; = 0.5 (ps)~^ (i.e., dot size / = 10 nm). 
As we can see, J'^^{uj) is comparable to JP^{uj) in that 
regime. But it shrinks much faster than Jp^{uj) as the 
dot size is increased and is negligible when the dot size 
I > 50 nm. Fig. 4 also shows the spectral functions 
JP^{uj) and J'^^{uj) at two different center-to-center dis- 
tances, with UJd = 0.05 (ps)~^ (i.e., d = 100 nm) and 
UJd — 0.02 (ps)"^ (i.e., d = 250 nm). The influence from 
the parameter d to both J^^{uj) and J'^^{uj) is small com- 
pared with that from the parameter /. It adds the spec- 
tral function an oscillation term, and the oscillation fre- 
quency (determined by ujd) is increased with d. All these 
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FIG. 4: Spectral functions of double quantum dot due to 
piezoelectric coupling (labelled by pz) and deformation po- 
tential coupling (labelled by df) with uJd = 0.02 and ujd ~ 
0.05 (ps)-^ (a) uji = 1 (ps)-^ (b) uJi = 0.5 (ps)-\ 



properties of spectral functions JP^{uj) and J'^^{uj) deter- 
mine the decoherence induced by piezoelectric coupling 
and deformation potential coupling to phonons, respec- 
tively. 

As one can see in Eq. (35) and Eq. (38), the deforma- 
tion potential coupling constant is real, while the piezo- 
electric coupling constant is imaginary, which means they 
do not interferei^^^ Thus the total spectral density is 



J{oj) = JP''{uj) + J'^^iuj) 



(40) 



B. Decoherence induced by acoustic phonons 

The decoherence of quantum system due to interacting 
with environment is a crucial point in quantum informa- 
tion. In a double quantum dot, scattering by phonons 
can cause considerable loss of coherence accompanied by 
dissipation when the tunneling electron flips back and 
forth between two dots. One of the advantages of our 
approach is that the decoherence rate in the this process 
is obtained explicitly. Thus one can analyze it clearly. 

Using the expressions of spectral density Jp^(uj) (Eq. 
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FIG. 5: Decoherence rates 7pz (solid line) and 7df (dash line) 
as functions of cutoff frequency loi when ujd ~ 0.02 (ps)~^. 
The tunneling rates in Fig. 5 (a) and Fig. 5 (b) are 0.05 (ps)~^ 
and 0.5 (ps)~^, respectively. 



(36)) and J'^^(cj) (Eq. (39)) above, the decoherence rates 
induced by piezoelectric and deformation potential cou- 
pling are written as 

2 2riTc LOd 



and 



7df = 277 gdfq^Tc^ {I - 



■sm )e 



, (42) 



respectively. Here, one should note the parameter 77 in 
Eq. (41) and that of the Eq. (42) are not the same, 
because they are calculated from Eq. (29) with different 
spectral functions (JP^(tj) and J'^^(a;), respectively). Ac- 
cording to Eq. (40), the total decoherence rate induced 
by acoustic phonons is 7 = 7pz + 7df . 

Fig. 5 presents the decoherence rates 7pz and 7df as 
functions of at two different tunneling rates = 
0.05 (ps)~^ and = 0.5 (ps)~^. Another parameter ujd 
is fixed as 0.02 (ps)""'^, which means the center-to-center 
distance between two dots is about 250 nm. As showed by 
Fig. 5(a), at small tunneling rate, the contribution to de- 
coherence rate arose from deformation potential coupling 




0.05 

Tc [(ps)-"] 



FIG. 6: Decoherence rate 7 as a function of tunneling rate Tc 
and distance between two dots d. The dot size is chosen to 
be 100 nm. 



is small compared with that from piezoelectric coupling, 
even in strong confinement regime (i.e., lj; ~ 1 (ps)~^). 
But when the tunnehng rate is large (Fig. 5(b)), the de- 
coherence rates arising from these two mechanisms are 
comparable, thus the contribution from the deformation 
potential coupling can not be neglected. Fig. 5(a) and 
(b) also show that the decoherence rates (both 7pz and 
7df) are suppressed when LUf, is decreased, which indicates 
that one should use large dot size to get little decoher- 
ence. However, large dot size means small characteris- 
tic energy spacing (on-site charging energy) wqi of single 
quantum dot. It is well known that, our two level Hami- 
tonian is valid to describe the double-dot system only if 
Tc, A:bT<Cwoi, where ks is Boltzmann constanti^ So, 
low temperature technique is needed to maintain good 
quantum properties of dots when the dot size is large, 
just as the experiment is performed. ^■^ 

In what follows, we choose a large dot size of 100 nm 
(approximate size for the dot in Ref. [13]), i.e., uji = 
0.05 (ps)^^. Since the tunneling barriers in experiment 
of Ref. [13] are made by depleting electrons with neg- 
ative gate voltage, their tunneling rates are flexibleiS 
In Fig. 6, we present the decoherence rate 7 (« 7pz at 
that dot size) as a function of tunneling rate and dis- 
tance between two dots (from 100 to 1500 nm). Some 
oscillations, coming from the sine term in the spectral 
density, are discerned from this 3-dimensional figure. We 
find the characteristic decoherence time T2 = I/7 spec- 
ulated from the figure is about 1 ns, which corresponds 
well with the value fitted from the experimental curved 
So the coupling to phonons is one of the main decoher- 
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ence mechanisms in such a double-dot system. It also 
shows that the decoherence rate increases with tunnel- 
ing rate when Tc < cui. But larger tunneling rate will 
suppress the decoherence drastically, even to zero when 
Tc w 0.1 (ps)-i (i.e., 2uji). This value of is stiU in 
the range of <C coai ~ 1.3 meV ~ 2 (ps)^^,^'^ in which 
our two level model holds. Thus, such kind of decoupling 
mechanism can be probably realized. 

V. CONCLUSION 

In conclusion, we have investigated the charge qubit 
dynamics in a semiconductor double quantum dot cou- 
pled to phonons at low temperature limit. Our approach 
is a perturbation theory after a unitary transformation. 
The dynamical tunneling current is obtained explicitly as 
a simple damped Rabi oscillation. The comparison with 
PER approach shows the advantages of our approach is 
that: it is not restricted by the form of spectral den- 
sity; it can be extended to strong coupling regime and 



works well for the whole range of tunneling rate T^, the 
long time behavior is also consistent with the symmet- 
ric double-dot system. Additionally, the simple deco- 
herence rate expression allows us to analyze the phonon 
induced decoherence clearly. We find that, in strong 
confinement regime of dot and large tunneling rate Tc 
(> 0.1 (ps)~^), the contribution to decoherence from de- 
formation potential coupling can not be neglected com- 
pared to that from piezoelectric coupling in GaAs mate- 
rial. The decoherence arose from both these two mech- 
anism will be suppressed when the dot size is increased. 
The decoupling with phonons will happen when the con- 
dition 2uji < Tc ujoi is met. 

Finally, we hope our predictions can be testified by 
experiment in the near future. 
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